Accepted to The Astrophysical Journal, April 6, 1999 
Preprint typeset using I^TgK style emulateapj v. 04/03/99 



LOW LUMINOSITY STATES OF THE BLACK HOLE CANDIDATE GX 339-4. I. 
ASCA AND SIMULTANEOUS RADIO/RXTE OBSERVATIONS 

JoRN Wilms', Michael A. Nowak^, James B. Dove^-\ Robert P. Fender"*, Tiziana Di Matteo^ * 

Accepted to The Astrophysical Journal, April 6, 1999 



Ph. 

< 

0^ 



> 
m 
(N 



ABSTRACT 

We discuss a series of observations of the black hole candidate GX 339-4 in low luminosity, spectrally hard 
states. We present spectral analysis of three separate archival Advanced Satellite for Cosmology and Astrophysics 
(ASCA) data sets and eight separate Rossi X-ray Timing Explorer (RXTE) data sets. Three of the RXTE obser- 
vations were strictly simultaneous with 843 MHz and 8.3-9.1 GHz radio observations. All of these observations 
have (3-9 keV) flux ^ 10"^ ergs s"' cm"^. The ASCA data show evidence for an w 6.4 keV Fe line with equiv- 
alent width w 40 e V, as well as evidence for a soft excess that is well-modeled by a power law plus a multicolor 
blackbody spectrum with peak temperature sa 1 50-200 eV. The RXTE data sets also show evidence of an Fe line 
with equivalent widths sa 20-140 eV. Reflection models show a hardening of the RXTE spectra with decreasing 
X-ray flux; however, these models do not exhibit evidence of a correlation between the photon index of the in- 
cident power law flux and the solid angle subtended by the reflector. 'Sphere-ndisk' Comptonization models and 
Advection Dominated Accretion Flow (ADAF) models also provide reasonable descriptions of the RXTE data. 
The former models yield coronal temperatures in the range 20-50 keV and optical depths of t « 3. The model fits 
to the X-ray data, however, do not simultaneously explain the observed radio properties. The most likely source 
of the radio flux is synchrotron emission from an extended outflow of size greater than 0(10^ GM/c^). 

Subject headings: accretion — black hole physics — Stars: binaries — X-rays: Stars 
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1. introduction 



The galactic black hole candidate (BHC) 
unique among persistent sources in that it shovl's 
ety of spectral states and transitions among thes( ; 
sumed order of increasing bolometric luminos ty, 
exhibits states with hard, po wer-law spectra ( 
vaisky et al. 1986; 'low state', iGrebenev et al 



GX 339-4 is 
a wide vari- 
: states. In pre- 
GX 339-4 
off state', Ilo- 
; a soft state 

with no e vidence of a power law tail ('high stite'; G rebenev 



19?1) 



et al. 1991); and a very bri ght, soft state with e) 



1951) 



law tail ('very high state'; Miyamoto et al 
are apparently times when the flux is high, but 
not as soft as the 'high' or 'very high state', 
der Klis (1997) refer to this as the 'intermediate 
note that there is some evidence of overlap 
Th e broad-band (GRANAT /SIGMA) hard state 
by Grebenev et al. (1991) apparently represent 
nous state than does the broad-band soft state 



the same instrument. [ fVliyamoto et al. (1995)| ha ^ 



possibility of hysteresis in galactic BHC state tr 
ilarly diverse sets of states h ave been observed 



sients such as Nova Muscae (Kitamoto et al. 1^ 92 



et al. 1994); however, GX 339-4 is closer to be 
source. 

Although there have been a number of olbservations of 



GX 339-4 in the near-infrared and op tical (Dox 
Motch et al. 1983|; iMotch et al. 1985|; Steiman 



1990; llmamura et al. 1990| ; [Cowley et al. 1991 



ameron et al 
including de 

tection of a 14.8 hr periodicity in the optical (pallanan et al 



tended power- 
There also 
he spectrum is 
Mendez & van 
state' . We also 
betvyeen the states, 
data presented 
i a more lumi- 
((ata taken with 
suggested the 
alnsitions.] Sim- 
in X-ray tran- 
; Miyamoto 
ng a persistent 



;ey et al. 1979 



1991 ), there is no convincing mass function for the system. In 
the o Jtical, the system is faint, variable (My ~ 16-20), and red- 
dene 1 (Ay = 3.5). The physical source of the optical emission 
is un mown. It has been hypothesized that it is entirely dom- 
inate 1 by the accretion disk, as the optical flux is apparently 
antic mela ted with the soft X-ray emission (S teiman-Cameron 
et al 1990; [imamura et al. 1990 ). These properties of the 
emission have made it difficult to obtain a good distance mea- 
surer lent, with estimates ranging from 1.3 kpc (P redehl et al. 
1991) 
4 kpc 
198^ 
by 



to 8 kpc 7 ^rindlay 1979 ), with many authors choosing 
( Doxsey et al. 1979|; C owley, Crampton & Hutchings 
. A careful study of these distance estimates is presented 
Z(|lziarski et al. (1998) who argue for a distance of « 4 kpc. 



G 



339-4 also has been detected in the radio (Sood & 



Cam )ell-Wilson 1994), and possib ly even has exhibited ex- 
tende d emission ( [Fender et al. 1997 ). Within the hard state, the 
spectrum is flat/inverted with a spectral index of a = 0. 1- 



ender et al. 1997; [Corbel et al. 1997 ), where the radio 



ct 



radio 
0.2 ( 
flux 
is 

nen 
trans 
comili, 
ley. 

Dilring the Rossi X-ray Timing Explorer (RXTE) Cycle 2 
obseiving phase (1996 December - 1998 February), we per- 
form ;d a series of eight RXTE observations of GX 339-4. The 
first hree observations were spaced a week apart from one 



lensity Sy ocv"". Furthermore, in this state the radio flux 
lated wi th the X-ray and gamma-ray fluxes (H annikain- 
al. 1998), but the radio flux disappears as GX 339-4 
ts to a higher X-ray flux/softer state (Fender 1998, priv. 
), which is comparable to the behavior of Cyg X-1 (P oo- 
ifender & Brocksopp 1998). 
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another from 1997 February 4 to 1997 February 18. These 
three observations were scheduled to be simultaneous with 8.3- 
9.1GHz radio observations that were conducted at the Aus- 
tralian Telescope Compact Array (ATCA ). The results of the 
radio observations have been reported by Corbel et al. (1997). 
Additionally, three 843 MHz observations performed at the Mo- 
lo ngolo Observatory Synthesis Telescope (MOST) and reported 
by Hannikainnen et al. (1998) are also simultaneous with these 
RXTE observations. 

This paper is structured as follows. We discuss the spectral 
analysis of archival Advanced Satellite for Cosmology and As- 
trophysics (ASCA) data in section ^. We look for evidence of 
Fe lines in the data and we characterize the soft (< 1 keV) X- 
ray data. In section H we present the RXTE data. We first dis- 
cuss the All Sky Monitor (ASM) data, and then we discuss the 
pointed observations. We perform spectral analysis much aki n 
to that which we considered for Cyg X-1 (Dove et al. 1998). 
Here, however, we consider Advection Dominated Accretion 
Fl ow (ADAF) models as w ell by using the models described 
by Di Matteo et al. (1998) We discuss the implications of the 
simultaneous radio data in section 0. In section ^ we discuss 
the implications of the X-ray observations for theoretical mod- 
els. We summarize our results in section |[ We present timing 
anal ysis of the RXTE data in a companion paper (N owak et al. 
1999, henceforth paper II). 

2. ARCHIVAL ASCA OBSERVATIONS 

The ASCA archives contain four observations of the 
GX 339-4 region. A log of the observations is given in in Ta- 
ble [l|. In appendix |^ we describe the methods that we used 
to extract, filter, and analyze these ASCA observations. To the 
best of our knowledge, an analysis of these observations has not 
been published previously, except fo r a power spectrum for on e 
of the observations (date not given, Dobrinskaya et al. 1997). 
The first of the observations (1993 September 16) did not de- 
tect the source, with the upper limit to the 3-9 keV flux being 
« 10"'^ ergs s"' cm^. As we will discuss further below, the in- 
ferred 3-9 keV fluxes for the remaining three observations (Ta- 
ble [l|) are lower by a factor of two to ten than the fluxes of the 
RXTE observations discussed in §|[ 

We chose to fit the ASCA data with a phenomenological 
model consisting of a multicolor blackbody spectrum plus a 
broken power law, considered with and without a narrow Gaus- 
sian line at « 6.4 keV. These fits are similar to those performed 
for ASCA observations of the hard state of Cyg X-1 (E bisawa 
et al. 1996), which shows evidence for a weak and narrow Fe 
line with equivalent width « 40 e V, as well as for a soft excess 
well-modeled as a multicolor blackbody with peak temperature 
w 150 eV. 

The fits with the phenomenological models yield rag- 
ing from 0.98 to 1 .4. The brightest data set showed the greatest 
evidence for structure beyond this simple model. A sample fit 
is shown in Figure Ju Note that the neutral hydrogen column 
was fixed to 6 x 10 cm"^. Allowed to freely vary, the neutral 
hydrogen column tended to float between 4 and 8 x 10^' cm"^, 
depending upon what combination of phenomenological mod- 
els was chosen, with minimal changes in the of the fits. As- 
sociated with these changes in best fit neutral hydrogen column 
were <; ±30% changes of the best fit peak temperature of the 
multicolor blackbody and even larger changes (factors of w 3) 
in the best fit normalization of the multicolor blackbody com- 
ponent. We should thus associate systematic error bars with 
these two parameters that are somewhat lai^ger than the statisti- 



cal error bars presented in Table ^ 

All fits improved with the addition of a narrow Gaussian line. 
In al l fits we fixed the line width to 0. 1 keV (see E bisawa et al. 
1996; who always found a < 0.1 keV in fits to ASCA data of 
Cyg X-1), and for the lowest flux data set we also fix the line 
energy to 6.4 keV. For the lowest flux data set the Ax^ = 5.5 for 
one addi tional parameter. By the F-test (B evington & Robin- 
son 1992), this is an improvement to the fit for one additional 
parameter at the 98% confidence level. The other two data sets 
show even more significant improvements to the x^- The best 
fit line equivalent widths ranged from « 30 to 60 e V. There is no 
compelling evidence for a strong flux dependence to the equiv- 
alent width of the line. 

The transition to the bright, soft state typically occurs at 3- 
9 keV luminosities > 10"^ ergs cm"^ s"' ; i.e., a factor of two to 
ten brighter than these ASCA observations. Thus these obser- 
vations offer useful tests of ADAF models, which are hypothe- 
sized to be mos t relevant to low-luminosity, hard stat e systems 
( [Narayan 199^ ; |Esin, McCUntock & Narayan 1997| ). ADAF 
models predict a detectable correlation between the tempera- 
ture of the soft excess, the strength of the Fe line, and the source 
luminosity. They hypothesize that the luminosity decay of BH 
transients is due, in part, to an increase of the radius at which the 
accretion flow transits from cold, geometrically thin, and radia- 
tively efficient to hot, geomet rically thick, and advective (E sin, 
McClintock & Narayan 1997). In some models, the transition 
radius can grow to as large as 0(10"* Rc), where Rq = GM/c^. 
(Such large t ransiti on radii are not a strict requirement of ADAF 
models; in §3.2.4 we show that somewhat smaller transition 
radii, « 200-400 Rc, are preferred for ADAF models of the 
RXTE data.) As discussed by Esin, McClintock & Narayan 
(1997), one then expects the peak temperature of the soft ex- 
cess to decrease below 150eV and the equivalent width of any 
Fe line to decrease to values less than « 30 eV. 

The best fit equivalent widths found for GX 339-4 are greater 
than can be accommodated in ADAF models with a large tran- 
sition radius, and they are also slightly larger th an pre dicted by 
the 'sphere-Hdisk' corona models de scribed in §3.2.3 (see also 
Dove et al. 1997; Dove et al. 1998). These latter models have 



a similar geometry to the ADAF models, and they often posit a 
coronal radius ^ 100 /?g- Likewise, we do not detect any large 
decreases in the best-fit disk temperatures with decreasing lu- 
minosity. Although it is dangerous to make a one-to-one cor- 
respondence between a phenomenological fit component and a 
physical component, these best-fit values are suggestive of, but 
not definitive proof of, temperatures hotter than can be accom- 
modated in models where cold, soft X-ray emitting material 
exists at very large radii. 

3. RXTE OBSERVATIONS 

3.1. The Monitoring Campaign 

To study the long-term behavior of GX 339-4, and to place 
our pointed observations within the context of the overall be- 
havior of the source, we used data from the All Sky Monitor 
(ASM) on RXTE. The ASM provides Ughtcurves in three en- 
ergy bands, 1.3-3.0keV, 3.0-5.0keV, and 5.0-12.2keV, typ- 
ically consisting of several 90s measurements per day (see 
Levine et al. 1996|; [Remillard & Levine 1997|; L ochner & 
Remillard 1997). In Figure ^ we present the ASM data of 
GX 339-4 up until Truncated Julian Date (TJD) « 1000 (1998 
July 6). We also indicate the dates of our RXTE observa- 
tions, as well as the dates of ATCA and MOST radio obser- 
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Table 1 

Log of the ASCA observations. 



Obs. 


Date 


Integration time 


SISO 


3-9 keV Flux 






(ksec) 


(cts s"') 


(10"^ ergs cm"^ s"') 


1 


1994 August 24 


15 


3.6 


0.11 


2 


1994 September 12 


17 


6.4 


0.19 


3 


1995 ScplemberS 


30 


17.5 


0.63 



Note. — All observations were taken in Bright 1-CCD mode. SISO: filtered SISO count rate. 



Table 2 

Parameters for a multicolor blackbody plus broken power law plus Gaussian line fits to ASCA data. 



Date 


Tin 


Adbb Ti 




Abpi 


El 


Al 


EW x^/dof 


Xred 




(keV) 


(xlO'') 


(keV) 


(xlO-^) 


(keV) 


(xlQ-*) 


(eV) 





1994 Aug. 24 




O 0+0.9 
^■^-1.6 


1 -7O+0.03 
^•'°-0.03 


o 4+0.4 


1 /:o+0.04 
i -0^-0.03 


^■^-0.1 








1500/1439 


1.04 


1994 Aug. 24 


n 14+""' 


^•^-1.6 


, -Tn+0.03 

i- 'O-0.03 


T Q+0.5 
^•-'-0.6 


1 <;/i+<'-03 
l-64_o.o3 


4 7+01 


6.4 


0.6_o.3 


34!?i 


1494/1438 


1.04 


1994 Sept. 12 




9 c+0.6 


1 1 +0.02 

i-°i-o.oi 


Q Q+0.2 

■^•°-0.2 




-7 9+0.1 
' -^-0.1 








1603/1621 


0.99 


1994 Sept. 12 




^•J-0.4 


1 Q 1+0.01 
i-oi-OOl 


o Q+0.1 

J-O-0.1 


1 «Q+0.01 
^•-'^-0.02 


-7 9+0.0 

' -^-0.1 


A o£+0.08 


1-0-0.2 


56!f 


1580/1619 


0.98 


1995 Sept. 08 


o.i9«:«? 


9 4-H).i 

^■^-0.3 


1 Qo-tO.02 
i-^J-OOl 


o ^+0.1 
J- '-0.1 


1 ^;n+o.oi 
ioO_o.o2 


25.6!^:3 








2597/1838 


1.41 


1995 Sept. 08 


0.19^S 


T 4-H).0 
^■^-0.1 


1 Ql+O.Ol 
•^•^^-0.00 


o -7+0.1 
J-'-OO 


i-Oj-OOl 


25.6!«:? 


"■-'i-0.06 


-J -j+0.6 


40!? 


2523/1836 


1.37 



Note. — Tin: peak multicolor blackbody temperature. Adbb^ multicolor blackbody normalization. Fi, Fz: broken power law photon indices. E\,: break energy. 
Abpi: Power law normalization (photons keV~' cm"^ s"' at 1 keV). Ey: hne energy. Ay: Line normalization (photons cm"^ s~' in the hne). EW: hne equivalent width. 
Uncertainties are at the 90% confidence level for one interesting parameter (Ax^ = 2.71). 




Fig. 1. — GX 339-4 ASCA observation of 1994 August 24 with energy bins rebinned by a factor of 25 for clarity. Model and associated 
residuals (data/model) are for the best fit multicolor blackbody plus broken power law without a Gaussian line component. For clarity, only the 
SISO (circles) and GIS2 (triangles) data are shown. 
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vations ( tl^ender et al. 1997|; [Corbel et al. 1997|; Hannikainnen 
et al. 1998). We discuss the long timescale variability of this 
lightcurve in paper II. 

Bas ed upon model fits to the observations of Grebenev et al. 



(1991) (low and high state), and [Vliyamoto et al. (1991) (very 
high state), we expect the different states of GX 339-4 to 
have ASM count rates as indicated in Table ||. The ASCA 
and RXTE observations discussed here are most characteristic 
of weak to average luminosity hard states. Confirmation that 
the eight RXTE observations taken between TJD 48 1 and 749 
do indeed represent a typical low/hard state comes from the 
broad band spectral analysis presented in §||, as well as from 
the timing analysis presented in paper II. The X-ray variabil- 
ity of these observations show root mean square variability of 
0(30%) and show a power spectrum (PSD) that, roughly, is 
flat below 0.1 Hz, cx /"' between 0.1-3 Hz, and oc above 
3 Hz. Time lags and coherence functions were also comparable 
to previously observed hard states of Cygnus X-1 (see paper II, 
and references therein). Further discussion and analyses of the 
timing data can be found in paper II. 

The transition to a higher flux level that occurs at TJD w 800 
appears to have a characteristic flux of a high state, but is 
not as soft in the 2-lOkeV bands as expected from the above 
cited high and very high states. This might be an example of 
what Mendez & van der Klis (1997) refer to as an 'interme- 



diate state' between hard and soft. No pointed RXTE obser- 
vations were taken during the transition, and four pointed ob- 
servations, which were not part of our monitoring campaign, 
occured shortly after the transition. Detailed confirmation of 
the spectral state suggested by the ASM data awaits analysis 
of these pointed observations. Note that the radio flux became 
quenc hed over the course of this transition to a higher ASM flux 
level (Hannikainnen et al. 1998). 



The variations observed in both the ASM lightcurve (prior 
to TJD w 800) and the pointed RXTE observations discussed 
below represent more than a factor five variation in observed 
flux. Comparable variations have been observed in the radio, 
and furthermore the radio lightcurves show evidence of a cor- 
relation with both the ASM a nd Burst and Transient S urvey 
Explorer (BATSE) lightcurves ( Hannikainnen et al. 1998| ). 



3.2. PC A and HEXTE Observations 

In this section we present the results from our analysis of the 
data from the two pointed instruments on RXTE: the Propor- 
tional Counter Array (PCA), and the High Energy X-ray Tim- 
ing Experiment (HEXTE). See appendix ^ for a description of 
the instruments and of the details of the data extraction and 
processing. A log of the RXTE pointed observations and the 
simultaneous radio observations is given in Table ^. 

As we show in appendix ^ there is a difference in the power- 
law slopes obtained from an analysis of spectra of the Crab with 
both instruments upon RXTE. In order to minimize the im- 
pact of this difference in the instrumental calibration onto the 
data analysis, we primarily analyze the data from both instru- 
ments individually and use the difference in model parameters 
between instruments as a gauge of the systematic uncertainties. 
We do perform some joint analysis of PCA and HEXTE data 
using various reflection models. In the following sections we 
discuss in detail the implications of the calibration uncertainty 
for our analysis. 

For our analysis of the RXTE broad band spectrum, we used 
several different spectral models consistent with the range of 
parameterizations currently used in the literature to describe 



the spectra of BHC. As in the ASCA analysis (§^, we fixed 
A^H = 6 X 10^' cm^. We first used the purely phenomenological 
exponentially cutoff power law and broken power law models 
as a broad character izatio n of the data. The results of this mod- 
eling are given in §3.2.1 and in Tables ^ and ^. We then ap- 
plied the three more physically motivated models that are cur- 
rently discussed in the literatu re: re flection of a power law off 
an (ionized) accretion disk (§3.2.2 an d Tabl e 'sphere and 
disk' corona Comptonization models (§ 3.2.3| and Table ||), and 
ADAF models (§3.2.4). The ADAF models are only applied 
to the unfolded data. Residuals for each of the best fits to the 
data from Observations 5 and 7 are shown in Figures^, ^ and^ 
(except for the ADAF models, where we present Observations 
1 and 5). We chose to present these former two observations 
because not only are they at extremes in terms of observed lu- 
minosity (Observation 5 is the faintest, and Observation 7 is 
the second brightest. Table Q), but also because they show de- 
tectable differences in their timing properties (paper II). 



3.2.1. Phenomenological Models 

Results from the purely phenomenological fits to the data, 
i.e., the broken power law and exponentially cutoff power law, 
are presented in Tables ^ and ||. We see that a broken power 
law plus a Gaussian models the PCA data very well. The low 
Xred (0.15-0.32) indicates that to some extent we may be fitting 
systematic features in the PCA response. The same may be 
true for the best fit parameters of the Fe line. The line widths 
(a K, 0.6 keV) and equivalent widths (« 130eV) are larger than 
for the ASCA observations, with the exception of Observation 5 
which has line parameters comparable to the ASCA observa- 
tions. As Observation 5 has the lowest count rate, it is more 
dominated by statistical errors and less dominated by system- 
atic errors than the other observations. Even ignoring the pos- 
sible systematic effects, however, we see that any observed line 
is narrower and weaker than is commonly observed in AGN. 

The « 3-10 keV spectral power-law slope is close to the 
'canonical value' of F = 1 .7; however, the PCA shows evidence 
for a hardening of this spectral slope above k, lOkeV. HEXTE 
data alone also show the high energy spectrum to be harder 
than the 3-lOkeV spectrum (Fig. |[ and Table |). Note that 
the difference between the PCA and HEXTE photon indices 
is greater than the discrepancy between the PCA and HEXTE 
best fit Crab photon indices (appendix and therefore it is 
unlikely to be a systematic effect. 

Such a hard HEXTE spectrum is consistent with previous 
observations by the Oriented Scintillation Spectrometer Exper- 
iment (O SSE) on board the Com pton Ga mma-ray Observatory 
(CGRO) (iGrabelskv et al. 19951 ; see also [Zdziarski et al. 1998| ). 
Grabelsky et al. (1995) found a slightly harder photon index of 
F = 0.88 and an exponential cutoff of = 68keV, somewhat 
lower than observed here. Note, however, that the OSSE ob- 
servations extended to k, 500 keV as opposed to the ~ 1 lOkeV 
of our HEXTE observation. Therefore, the HEXTE data for 
GX 339-4 do not strongly constrain the exponential rollover, 
and slightly harder power laws with lower exponential cutoffs 
are permitted. 



3.2.2. Reflection Models 

A spectral hardening above k, 7keV is the expected sig- 
nature of reflection of a hard continuum off of cold material 



(Magdziarz & Zdziai'ski 1995 ). Ueda, Ebisawa & Done (1994) 



applied reflection models to Ginga data of GX 339-4 and found 
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200 400 600 800 1000 200 400 600 800 1000 

JD - 2,450,000.5 JD - 2,450,000.5 

Fig. 2.— Left: RXTE All Sky Monitor data for GX 339-4 (5 day averages in the 1.3-12.2keV band) vs. Truncated Julian Date (TJD) = 
Julian Date (JD) -2450000.5. Dashes indicate dates of our RXTE pointed observations, diamonds indicate dates of MOST radio observations 
(Hannikainnen et al. 1998), and triangles indicate dates of ATCA radio observations (Fender et al. 1997, Corbel et al. 1998). Right: ASM colors 
for GX 339-4 vs TJD. Colors shown are: (1.3-3 keV)/(3-5 keV), (3-5 keV)/(5-12 keV)+5, and (1.3-3 keV)/(5-12.2 keV)+10 (the latter two ratios 
have been offset for clarity). 



Table 3 

Approximate expected ASM colors for the different states of GX 339-4. 



State 


count rate 


1.3-3 keV/ 


1.3-3 keV/ 




cts s"' 


3-5 keV 


5-12.2keV 


low 


7 


2 


2 


high 


15 


4 


30 


very high 


60 


4 


15 



Table 4 

Log of RXTE AND RADIO OBSERVATIONS. 



Obs 


Date 


TJD 


T 


Rate 


3-9 keV 


9-30 keV 


30-100 keV 


ATCA 


MOST 


a 








(ksec) 


(cps) 




(10"' erg cm"' s"') 




(mJy) 


(mJy) 




01 


1997 Feb. 02 


481 


11 


830 


1.07 


1.68 


2.65 


9.1 ±0.2 


7.0±0.7 


0.11 


02 


1997 Feb. 10 


489 


10 


730 


0.94 


1.50 


2.41 


8.2±0.2 


6.3±0.7 


0.11 


03 


1997 Feb. 17 


496 


8 


700 


0.90 


1.43 


2.35 


8.7±0.2 


6.1 ±0.7 


0.15 


04 


1997 Apr. 29 


567 


10 


470 


0.60 


0.97 


1.55 








05 


1997 Jul. 07 


636 


10 


200 


0.25 


0.43 


0.75 








06 


1997 Aug. 23 


683 


11 


650 


0.74 


1.18 


1.98 








07 


1997 Sep. 19 


710 


10 


730 


0.96 


1.48 


2.36 








08 


1997 Oct. 28 


749 


10 


480 


0.63 


1.01 


1.68 









Note. — We list: T, the duration of the RXTE observations; the average PCA count rate; the average (3-9 keV), (9-30 keV), and (30-100 keV) energy fluxes 



(all normalized to the PCA calibration); the flux density of fhp R ^-Q 1 hHt atc 



: Alogi^/AlogSy, the spectral index of the radio observations (Fender et al. 1997 
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Corbel et al. 199' 
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fliiT Hmicity nf the MHz MOST observations; and 



Hannikainnen et al. \99i ). 
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Fig. 3. — Spectral modeling of the PCA data from Observation 5 (left) and 7 (right). Residues are shown as the contribution to x- a) Count rate 
spectrum and the best-fit broken power-law with Gaussian line, b) Contribution to x from the brolcen power-law with Gaussian line, c) Contribu- 
tion to X from the ionized reflection model (pexriv) with Gaussian line, d) Contribution to x from the sphere-i-disk model with additional Gaussian 
line. 



Table 5 

Parameters for Gaussian line plus broken power law models (PCA only). 
Model: exp(-A'HiTa)(gauss(£;ii„e,iT,Aime)+bknpower(ri,r2,£breiik,Abim)) 
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Fig. 4. — Spectral modeling of the HEXTE data from observation 5 (left) and 7 (right). Residues are shown as the contribution to x- a) Count 
rate spectrum and the best-fit power-law with exponential cutoff, b) Contribution to x from the power-law with exponential cutoff, c) Contribution 
to X from the sphere+disk model. 



Table 6 

Parameters for exponentially cutoff power law models (HEXTE only). 



Model: const. (ApL^'^exp {-E/Eout}) 
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Fig. 5. — Spectral modeling of the joint PCA and HEXTE data for Observation 5 and 7. Residues are shown as the contribution to X- a) Count 
rate spectrum and the best fit broken power-law with Gaussian line, b) Contribution to x from the best fit broken power-law with a Gaussian line 
(parameters not given in text), c) Contribution to x from the best fit ionized reflector model with a Gaussian line. 
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strong evidence of reflection, whereas 3rabelsky et al. (1995) 



found no evidence of reflection in OSSE data. Z dziarski et al. 
(1998) jointly fit these simultaneously observed data sets and 
find that reflection models, albeit with a large Fe abundance, 
provide a ver y good description of the data. We have applied 
the models of [Vlagdziarz & Zdziarski (1995) , as implemented 
in XSPEC ipexrav, pexriv), to the GX 339-4 data. These mod- 
els consider an exponentially cutoff power law reflected off of 
neutral (pexrav) or partially ionized {pexriv) cold material. 

In Table ^ we show the fit results for reflection off of partially 
ionized ma terial similar to the models presented by Z dziarski 
et al. (1998). Just as in [Zdziarski et al. (1998)| , we include a disk 
component where we fix the inner disk temperature to 250 eV. 
As PCA does not usefully constrain models below 3 keV, the 
disk component is not strongly constrained; typically the 
values were higher by 5-20 without this component. We also 
fix the reflector inclination angle at 45°, fix the disk temper- 
ature at Tdisk = 10^ K, freeze the abundances at solar, but let 
the Fe abundance be a free parameter. In all our fits we found 
that the Gaussian line width, a, would tend to drift toward 0, 
so we fixed a = 0.1 keV. For the combined PCA and HEXTE 
data, we also fixed the Gaussian line energy to 6.4 keV. For 
fits to PCA data alone and joint PCA/HEXTE data, the expo- 
nential cutoff energy, Em^, would drift towards very large en- 
ergy (^ 1000 keV). We therefore only consid ered pure power 
laws without cutoffs. Zdziarski et al. (1998) have argued that 
the high energy cutoff is sharper than exponential, which one 
would not expect to be strongly constrained by the combined 
PCA/HEXTE data. 

A s for the Ginga data of GX 339-4 (U eda, Ebisawa & Done 
1994), the PCA data alone were extremely well described by 
reflection models. Again, however, the extremely low x^d (^^ 
low as 0.28) makes us caution that these models might partly 
be fitting systematics in the PCA response. PCA and Ginga 
are also very similar instruments in terms of design, and so to 
some extent they should exhibit similar systematic effects (as 
discussed in appendix ^ the internal consistency of the PCA 
calibration is now as good as or better than that for the Ginga 
calibration.). Note that the best fit Fe line equivalent widths 
here are significantly smaller than those found with the purely 
phenomenological models discussed in § [3.2. l| . 

The fits for the HEXT E data alone (not shown) were simi- 
lar to flie OSSE results of prabelsky et al. (1995)[ Namely, if 
one allows for an exponential cutoff (typically w lOOkeV) to 
the power law, the best fit reflection fraction becomes / < 0.01 . 
Such a small reflection fraction is not surprising considering 
how well a pure exponentially cutoff power law fits the HEXTE 
data (Table^. If one does not allow an exponential cutoff, the 
reflection fraction becomes / 3. Such a fit is trying to mimic 
a hard power law with a high energy cutoff. 

A joint analysis of the PCA and the HEXTE data should be 
similar to a joint analysis of the Ginga and OSSE data. Indeed, 
suc h an analysis yields re sults comparable to those presented 
by Zdziarski et al. (1998) if we constrain the photon index of 
the incident power law to be the same for both the PCA and 
HEXTE data. Notable for the results of such fits (not presented) 
is the fairly large overabunda nce of Fe (Ape = 3.2-5.2). Simi- 
larly, Zdziarski et al. (1998) find a large Ape = 2.5-3.0 except 
for a short data set, more likely dominated by statistical errors 
rather than systematic errors, where they find Ape = 1.6-2.0. 



For our joint PCA/HEXTE data, the best-fit reflection fractions 
were approximately 20% larger than the best-fit reflection frac- 
tion for PCA data alone. Such an increase in reflection frac- 
tion in general will reproduce the spectral hardening seen in the 
HEXTE energy bands. Increasing the average best fit Fe abun- 
dance from (Ape) = 1.6 (PCA data alone) to (Ape) =4.0 (joint 
PCA/HEXTE data) also leads to an increased spectral harden- 
ing above w 7 keV, while leaving the spectrum below sa 7 keV 
relatively unchanged. That is, such a fit helps to reproduce the 
spectral break at w lOkeV. 

For the joint PCA/HEXTE analysis, there is clearly a worry 
that these results are influenced by the systematic differences 
between the PCA and the HEXTE responses. We therefore per- 
formed reflection model fits where we constrained all fit param- 
eters to be the same for the PCA and the HEXTE data except for 
the incident power law photon index, which we allowed to vary 
between the two instruments'. Such models provided reason- 
ably good fits to the data, with Xred ranging from 0.67 to 1.09. 
The difference between the PCA and the HEXTE best fit pho- 
ton indices ranged from 0.05 to 0.13, with an average value of 
0.08. This is consistent with the systematic difference between 
the best-fit photon indeces for the Crab spectrum. For these 
models we find (Ape) = 2.3, which is more consistent with the 
results for fits to the PCA data only, and is sl ightly smaller than 
the results found by Zdziarski et al. (1998)| Note that we also 
find smal ler values of the ioniz ation parameter, ^, than were 
found by [Zdziarski et al. (1998)| . 

3.2.3. Corona Models 

We con sidered 'spheren-disk' Comptonization models (D ove 
et al. 1997) of the GX 339-4 observations. We have previously 
applied these mode ls to an RXTE observation of Cygnus X-1 
(Dove et al. 1998). The models consist of a central, spheri- 
cal corona surrounded by a geometrically thin, flat disk. Seed 
photons for Comptonization come from the disk, which has a 
radial temperature distribution ^7disk(^) c>c R~^l'^ and a temper- 
ature of 150 eV at the inner edge of the disk. Hard flux from 
the corona further leads to reflection features from the disk or 
to soft photons due to thermalization of the hard radiation. The 
(non-uniform) temperature and pair balance within the corona 
are self-c onsistently calculated from the radiation field (D ove 
et al. 1997). 

we parameterize our 



As described by Dove et al. (1997) 



models by the coronal compactness 

(Tt Lc 



4 = 



meC^ Rc 



(1) 



where crj is the Thomson cross section, me is the electron 
mass, Lc is the luminosity of the corona, and Rq is the ra- 
dius of the corona. Likewise, we define a disk compactness, 
= {\- ft:){aj / m^c^)PG / Rc, where Pq is the total rate of grav- 
itational energy dissipated in the system, and /c is the fraction 
dissipated in the corona. In calculating the numerical models, 
we set f d = 1 ■ Models with other values of 1^ yield the same 
ranges of self-consistent c oronal temperature s and opacities. In 
general = 4/(4 + 4) ( [Dove et al. 1997[ ). Based upon the 
'sphereH-disk' geometry, a fraction f 0. 32 of the coronal flux 
is absorbed by the disk ( Dove et al. 1997 ). The models are fur- 
ther parameterized by an initial electron coronal optical depth, 

'The photon index was constrained to be the same for HEXTE Cluster A and B. The necessary different normalizations between the PCA and the HEXTE models 
were subsumed into the constants multiplying the HEXTE models. As HEXTE requires a harder power law. these constants were now 0.42-0.57, as opposed to 
~ 0.7. Furthermore, the constants showed larger uncertainties, as the uncertainty of the HEXTE photon index now couples strongly to the value of the constants. 
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Tc (approximately equal to the total optical depth, as pair pro- 
duction is negligible for the parameters of interest to us), and 
a normalization constant Akot- From the best-fit compactness 
and optical depth, the average coronal temperature can be cal- 
culated a posteriori. 

Attempts to fit these models to the joint PCA/HEXTE data 
failed. Typical Xred values, even allowing for the inclusion 
of an extra Gaussian line component, were > 1.3. These fits 
showed a clear tendency for a hardening in the HEXTE band, 
and therefore we considered them to be influenced by the cross- 
calibration uncertainties between the PCA and the HEXTE in- 
struments (Note that our previous fits of Cyg X- 1 used an earlier 
version of the PCA response where we applied 1.5% system- 
atic uncertainties across the entire PCA band; these fits yielded 
Xred ~ 1-6' without considering an additional Gaussian compo- 
nent). We therefore considered 'sphere-ndisk' models fit to the 
PCA and the HEXTE data separately. In Table ^, we present the 
best-fit parameters for these models applied to our GX 339-4 
data. 

Although our numerical 'sphere-ndisk' models do include 
reflection and a fluorescent Fe line (typical equivalent width 
« 25 eV) from the disk, the PCA data showed residuals in the 
5-7 keV band, sim ilar as in in our fits to RXTE data of Cyg X-1 
( Dove et al. 1998| ). We included an additional Gaussian com- 
ponent to our fits. The equivalent widths of the additional lines 
were « 150 eV, and they appeared to be broad {a « 0.8 keV). 
This additional line may be attributable partly to uncertainties 
in the PCA response. For these fits, as well as for the reflection 
model fits, lines with energies significantly less than 6.4 keV 
can be fit, and this is likely a systematic effect. Part of the dis- 
crepancy between the data and the model, however, is signifi- 
cant. As we have discussed for our fits to the RXTE Cyg X-1 
data, there are several possible physical interpretations for the 
additional required equivalent width: There may be an over- 
lap between the disk and the sphere (our models invoke a sharp 
transition), the disk may be flared (we model a flat disk), the 
disk may have non-solar abundances, or alternatively one might 
invo ke a 'patchy disk' embedded in the corona (Z dziarski et al. 
1998). The best-fit reflection fractions of / « 0.4-0.5 found 
above are further indication that our models may require an ad- 
ditional source of reflected flux. 

Allowing an additional Gaussian line component, the fits to 
the PCA data yield extremely low x^d' which could be indi- 
cating that we are partly fitting systematic features in the PCA 
response. Note also that the PCA data fits yielded consistently 
larger optical depths and consistently lower compactness pa- 
rameters than the HEXTE data fits. The latter was more sig- 
nificant, and is again indicative of the HEXTE response being 
harder than the PCA response. Both instruments yielded op- 
tical depths Tc « 3^; however, due to the discrepancy in the 
best-fit spectral slopes between the PCA and the HEXTE bands, 
the best-fit average coronal temperatures range from 2 1-30 keV 
(PCA) to 34^5 keV (HEXTE). 

3.2.4. ADAF Models 

The basic picture of mass accretion via an ADAF in the c on- 
text of galactic BHC was introduced by Ichimaru (1977) and 
has been elab orated upon in a series of papers by Narayan and 
collaborators ( Narayan, Kato & Honma 1997|; E sin, McClin- 
tock & Narayan 1997). The accretion flow is divided into two 
distinct zones: the inner part is modeled as a hot, optically thin 
ADAF similar in some respects to the spherical corona dis- 
cussed above, while the outer part consists of a standai^d op- 



tically thick, geometrically thin disk. The transition radius be- 
tween the two zones, rt, = R,r/Rc, is one of the model parame- 
ters. We compute the ADAF spectrum acc ording to the proce- 
dure described by pi Matteo et al. (1998)t The electrons in an 
ADAF cool via three processes: bremsstrahlung, synchrotron 
radiation, and inverse Compton scattering. In addition we add 
the emission from a thin disk — calculated as a standard mul- 
ticolor blackbody — and include the Compton reflection com- 
ponent due to the scattering of high energy photons incident on 
the disk. 

In the ADAF models discussed here, we fix the black hole 
mass to be m = M /Mq = 6, assume the magnetic field to be in 
equipartition with thermal pressure {(3 = 0.5), and set the stan- 
dard Shak ura-Sunyaev viscosity parameter (S hakura & Sun- 
yaev 1973) to be ass = 0.3. We normalize the accretion rate 
to m = Mc^/LEdd, where Lgdd is the Eddington luminosity of 
the source. The hard state corresponds to mass accretion rates 
m < ;7ici-it = 10"^, where nicnt is the critical value above which 
an ADAF no longer exists. As m increases towards merit, the 
scattering optical depth of the ADAF goes up which causes the 
spectrum to become harder and smoother. Most of the flux from 
the ADAF plus disk configuration is emitted around lOOkeV 
and the spectrum falls off at higher energies. 

The model spectrum changes mainly as a function of r^- and 
m. The various spectral states correspond to different values of 
these parameters. For example, E sin, McClintock & Narayan 
(1997) attempt to explain the initial transition from soft to hard 
seen in the decay of Nova Muscae by a large change in r^ (from 
rti « 10 to rti « 10"*), followed by an exponential decay in m for 
the subsequent evolution of this transient system. The ASCA 
data of GX 339-4 discussed in §^ imply that comparably large 
changes in r^i are not relevant to those observations. Here, how- 
ever, unfolded RXTE data from Observation 1 and Observa- 
tion 5, the brightest and faintest observation respectively, can 
be described by ADAF models with r^ = 200, m = 0.08 and 
rti = 400, m = 0.05, respectively. These model spectra and 
RXTE spectra for Observations 1 and 5 unfolded with a cut- 
off broken power law plus Gaussian line are shown in Fig. |[ 
In these ADAF models, the observed spectral and luminosity 
changes of GX 339-4 are predominantly driven by changes in 
the transition radius; the implied accretion rate change is sub- 
stantially smaller than the factor of 5 for the observed luminos- 
ity change. 

For rti ^ 200, which provides a rough description of Obser- 
vation 1, the disk blackbody emission peaks in the far UV/soft 
X-rays and dominates over the s ynchrotron emission, which 
peaks in the optical/UV (see also Mziarski et al. 1998). The 
peak synchrotron emission frequency scales as oc rrT^^^rh^^^r^^^ 
and peaks in the range ly = 10"-10'^ Hz for supermassive black 
holes and i/ = lO'^-lO'^ Hz for galactic black holes. The spec- 
trum below the peak is approximately S^o^v^- The synchrotron 
emission can contribute significantly to the radio emission o f 
super-massive black holes (although see Di Matteo et al. 1998 ); 



however, the predicted radio flux of GX 339-4 is ten orders of 
magnitude below the observed 7 mJy flux at 843 MHz. Thus, 
there must be an extended source of radio emission, which we 
further discuss in the next section. 

4. SIMULTANEOUS RADIO OBSERVATIONS 

The first three of our RXTE observations were simultane- 
ous with 843 MHz observations taken with the Molongolo Ob- 
servatory Synthesis Telescope (MOST), and with 8.3-9.1 GHz 
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Table 8 

Models of 'sphere+disk' Comptonization plus a Gaussian line fit to PCA data only and HEXTE data only. 
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Note. — The three fit parameters of the Comptonization model are the compactness of the corona, £c, the coronal optical depth, Tc, and a normalization constant, 
Ajco,. The Gaussian line is parameterized as in the previous tables. From the best fit parameters, the equivalent width of the Une, EW, and the density averaged coronal 
temperature, kTc, are derived. 




Fig. 6. — Advection Dominated Accretion Flow models for the unfolded RXTE data from Observation 1 (solid line) and Observation 5 (dash-dot 
line). Parameters consistent with the unfolded data are described in the text. A source distance of 4 kpc was assumed. 
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observations taken at the Australian Telescope Compact Array 
(ATCA). Extensive disc ussion of the MOST and AT CA obser- 
vations can be found in [Hannikainnen et al. (1998) and Corbel 
et al. (1997), respectively (see also Table 

An estimate of the minimum size of the radio emitting re- 
gion can be obtained by noting that observationally the bright- 
ness temperatures of radio sources usually are not larger than 
10'^ K, else the electrons will suffer catastrophic inverse Comp- 
ton losses. The brightness temperature of a uniformly bright 
spherical source is given by {cD / dvf'S^ /lirk, where d is the 
diameter of the source, D is its distance, v is the observed radio 
frequency. Si, is the observed radio flux density, c is the speed 
of light, and k is the Boltzmann constant. Taking the 7 mJy 
observed at 843 Mz by MOST and the fact that GX 339-4 is 
unresolved, we derive 



£/>4x lO'^cm 



/ D 



\4kpc 



!3x10*'7;g' ^ 



V4kpc J VlOM, 



M 



(2) 



which is orders of magnitude larger than the inferred size of 
the X-ray emitting regio n, even for models that posit extremely 
extended coronae (e.g., p3sin. McClintock & Narayan 1997 ; 



Kazanas, Hua & Titarchuk 1997 ) 

This size scale is oc so that emission at 8.6 GHz could 
arise in a region an order of magnitude smaller than that respon- 
sible for the emission at 843 MHz. Indeed, the flat spectrum 
emission is likely to arise in a conical jet with a radia lly de- 
creasing optical depth (e.g., Hjellming & Johnston 1988 ). Thus 
the outflow likely has an extent of 0(10' GM/c'^) or greater 
Similar estimates for source size have been made for the other 
persistent b lack hole candidate and Z -source neutron star X-ray 



binaries by Fender & Hendry (1999) 



Assuming a radio spectral index of a = 0.1 and a sharp cut- 
off at 10/im (a reasonable upper limit for where the radio flux 
becomes optically thin, and typical of where ADAF models be- 
come optically thin in the radio), the radio flux is approximately 
0.1% of the 3-100 keV X-ray flux. The correlation between the 
X-ray and radio fluxes found by Hannikainnen et al. (1998) , 
co mparable to the X-ray/radio correla tion observed in Cyg X- 
1 ( l^ooley. Fender & Brocksopp 1998 ), suggests that there is a 
coupling between the inner accretion disk and the extended out- 
flow on timescales of 7 days or less. Matter leaving the corona 
at the escape velocity (0.25c at 3QRg) and thereafter decelerat- 
ing under the influence of gravity would take roughly 7 days to 
travel a distance of 10^/?g- As 7 days is the upper bound to the 
correlation timescale, the radio emitting outflow must leave at 
slightly greater than escape velocity, or there must be at least 
some amount of acceleration of the outflow. 

Although the radio observations are strictly simultaneous 
with our first three RXTE observations, GX 339-4 exhibits less 
than 1 % rms variability over the shortest timescales for which 
a reasonable radio flux estimate can be made (> 10 minutes). 
Thus there are no strong features to correlate between the radio 
and X-ray bands. 

5. DISCUSSION 

Coronal Size and Luminosity Variation: 

The relationship between the inferred size of the corona and 
the magnitude of the observed flux depends upon w hich s pec- 
tral model we are considering. As discussed in §3.2.4, for 



ADAF models one can associate lower fluxes with increased 
coronal radii. A larger coronal radius implies a lower efficiency 



and hence a decreased observed flux, even for constant accre- 
tion rates. Paper II shows that the characteristic power spectral 
density (PSD) timescale for GX 339-4 decreases for the low- 
est observed flux (Observation 5). If one associates the PSD 
timescale with characteristic disk timescales, this could be in 
agreement with an increased coronal radius. However, in pa- 
per II we also show that the time lags between hard and soft 
X-ray variability decreases with decreasing flux, which seems 
counter to a positive correlation between flux and coronal size. 

The 'sphereH-disk' coronal models make no assumptions 
about the radiative efficiency of the accretion. The flux can 
be either positively or negatively correlated with coronal ra- 
dius, depending upon the variations of the coronal compactness, 
4, and the temperature, T^, at t he inner edge of t he accretion 
disk that surrounds the corona (Dove et al. 1997). Note that 



the 'sphereH-disk' models used in this work, contrary to many 
ADAF models, do not consider synchrotron photons as a source 
of seed photons for Comptonization. 

Using the definitions of 4, ^i, and / given in §3.2.3, energy 
balance in the 'sphere-ndisk' system determines the coronal ra- 
dius, to within factors of order unity, to be given by 



Rc ~ 160 



W4 
4+4 



1/2 



D 



4kpc 



V 150 eV 



6M(r 

M 



10 '^ergs cm ^ s 



1/2 



Rg 



(3) 



where M is the mass of the compact object, D is the distance to 
the source, and Ftot is the bolometric flux of the source. If 7d, 
/, li, and 4 were held fixed, then the coronal radius would be 
positively correlated with flux. Whereas this might pose some 
problems for understanding the flux dependence of the charac- 
teristic timescales observed in the PSD, this would agree with 
the flux dependence of the X-ray variability time lags (paper II). 
However, as the RXTE bandpass does not usefully extend be- 
low w 3 keV, we do not have a good handle on the flux depen- 
dence of Ta. If 7d oc fjgj with /? > 1/4, then increasing flux 
could imply decreasing coronal radius. 

Correlations Among Spectral Para meters: 

ueaa, HDisawa & uone (iyy4y claimed that reflection mod- 
els of GX 339-4 exhibited a correlation between photon in- 
dex, r, and reflec tion fraction, f, w ith softer spectra implying 
greater reflection. Zdziarski (1998) has claimed that this corre- 
lation extends to reflection models of Seyfert 1 galaxies as well. 
Such a correlation is not unreasonable to expect. For example, 
if we allow the corona a nd disk to overlap to some ext ent in 
the 'sphereH-disk' model ( ^outanen, Krolik & Ryde 1997 ), then 
we expect the increase in the flux of seed photons to cool the 
corona and lead to a softer spectral index. Likewise, the cover- 
ing fraction of the disk would be increased, in agreement with 
the suggested correlation. In Figure |^ we plot / vs. F for our 
reflection model fits to GX 339-4. Contrary to the claims of 



Ueda, Ebisawa & Done (1994) and [Zdziarski (1998) , however, 
there is no strong evidence for a correlation. Fitting the reflec- 
tion fraction with a function linear in F, as opposed to fitting 
with the mean value of /, improves the of the fits by 0.2, 
which is not significant. Fitting with the mean gives Xred - 

We do note two possible trends from the reflection model 
fits. First, as has been noted for other h ard state galactic black 
hole candidates (Fanaka & Lewin 1995 and references therein). 



there may be a correlation between flux and photon index with 
lower flux implying a harder source. (The significance of the 
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Fig. 7. — Left: Reflection fraction vs. photon index, F, for models fit to PCA data only (squares), and models fit to PCA plus HEXTE data 
that allowed the PCA and HEXTE photon indeces and normalizations to be different (diamonds; HEXTE photon index shown). Middle: Photon 
index, F, vs. observed 3-9 keV flux for the same reflection models as on the left. Also shown, without error bars, is the best fit compactness, (.c, 
for 'sphere+disk' coronal models fit to HEXTE data (small triangles). Right: Disk ionization parameter, ^, in units of ergs cm s"' for the same 
reflection models as on the left. 



correlation is driven by Observation 5, the faintest and hard- 
est of the observations. However, a similar correlation is also 
present in color-intensity diagrams.) Such a correlation is con- 
sistent with the expectations of ADAF models where the ra- 
dius increases with decreasing flux (Figure Again, the 
'sphereH-disk' corona models do not predict a clear trend with- 
out knowing the flux dependences of other parameters such as 

Second, the ionization parameter, ^, is positively correlated 
with flux. S uch a correlation was noted by Z ycki, Done & 
Smith (1998) for Ginga observations of Nova Muscae. It is 
not unreasonable to expect the disk to become increasingly ion- 
ized with increasing flux. We again caution that it is danger- 
ous to make one-to-one correlations between a model fit pa- 
rameter and a true physical parameter Furthermore, the sig- 
nificance of the correlation is again almost entirely determined 
by Observation 5, the faintest observation, which has ^ ~ 0. 
However, if we take the flux dependence of ^ as being real 
and interpret it physically, it provides some constraints on the 
flux dependence of the coronal radius. The ionization param- 
eter is cx Fxot/ipR^), where p is the density of the disk. For 
a gas pressure-dominated S hakura-Sunyaev a-disk, p oc 
(Shakura & Sunyaev 1973). In order for ^ to be roughly lin- 



ear in flux (the actual dependence is not strongly constrained 
by the data), we require Ftot K R^'^^ . (For the 'sphere-ndisk' 
models this would further require 7d PS 7?"'/-', depending upon 
the flux dependences of l^, 4, g, etc.) Taken physically and 
in the context of a gas pressure dominated Shakura-Sunyaev a- 
model, the flux dependence of ^ implies that the coronal radius 
increases with increasing flux. 

6. SUMMARY 

We have presented a series of observations of the black hole 
candidate GX 339-4 in low luminosity, spectrally hard states. 
These observations consisted of three separate archival ASCA 
and eight separate RXTE data sets. All of these observations 
exhibited (3-9 keV) flux < 10"^ ergs s"' cm"^, and the observed 
fluxes spanned roughly a factor of 5 in range for both the ASCA 
and RXTE data sets. Subject to uncertainties in the cross cali- 
bration between ASCA and RXTE, the faintest ASCA observa- 
tion was approximately a factor of two fainter than the faintest 
RXTE observation. 



All of these observations showed evidence for an w 6.4 keV 
Fe line with equivalent widths in the range of k, 20-140eV. 
The ASCA observations further showed evidence for a soft 
excess that was well-modeled by a power law plus a multi- 
color blackbody spectrum with peak temperatures in the range 
« 150-200eV. Both of these factors considered together argue 
against 'sphere-ndisk' or ADAF type-geometry coronae with 
extremely large coro nal radii of 0(10"* Rg) (e.g., E sin, McClin- 
tock & Nai-ayan 1997). 

The RXTE data sets were well-fit by 'sphere-ndisk' Comp- 
tonization models with coronal temperatures in the range 20- 
50 keV and optical depths in the range of r w 3. These fits were 
similar to our previous fits to RXTE data of Cyg X- 1 . Advec- 
tion Dominated Accretion Flow models, which posit a similar 
geometry, also provided reasonable descriptions of the unfolded 
RXTE data. The 'sphere-ndisk' and ADAF models were not 
able, however, to also model the observed radio fluxes. Thus, 
a static corona seems to be ruled out by the observations. The 
ADAF models can imply that the coronal radius increases with 
decreasing flux. The 'sphere-ndisk' corona models do not make 
a specific prediction for the dependence of the coronal radius on 
the flux; however, they can be consistent with a positive correla- 
tion between coronal radius and flux. As described in paper II, 
a positive correlation between flux and coronal radius is consis- 
tent with the observed flux dependence of the time lags between 
hard and soft X-ray variability. 

We also considered 'reflection models' of the RXTE data. 
These models showed evidence of a hardening of the RXTE 
spectra with decreasing X-ray flux. They further showed evi- 
dence of a positive correlation between the best-fit ionization 
parameter, ^, and the observed flux. Especially the latter of 
these correlations, however, was dominated by the model fits of 
the faintest observation. The reflection models did not exhibit 
any evidence of a correlation between the photon index of the 
incident power law flux and the solid angle subtended by the 
reflector 

Three of the RXTE observations were strictly simultaneous 
with 843 MHz and 8.3-9.1 GHz radio observations. The most 
likely source of the radio flux is synchrotron emission from an 
extended outflow with a size of 0(10^ GM/c^). The correla- 
tion bet ween radio and X-ray emi ssion on timescales of 7 days 
or less ( Hannikainnen et al. 1998) implies a strong coupling of 



Wilms et al. 



15 



the inner disk accretion flow with this spatially extended out- 
flow as is expecte d by recent theoretical arguments (B landford 
& Begelman 1998). Further simultaneous radio/X-ray observa- 
tions, preferably with the addition of IR/optical monitoring to 
constrain the location of the synchrotron break and with the ad- 
dition of soft X-ray monitoring to constrain the accretion disk 
parameters, are required to test such models in detail. 
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APPENDIX 
ASCA DATA EXTRACTION 

We extracted data from the two solid state detectors (SISO, SISl) a nd the two gas d etectors (GIS2, GIS3) onboard ASCA by using 



the standard ftools as described in the ASCA Data Reduction Guide ( [Day et al. 1998[ ). The data extraction radius was limited by the 
fact that all the observations were in 1-CCD mode and that the source was typically placed close to the chip edge. We chose circular 
extraction regions with radii of w 4 arcmin for SISO, « 3 arcmin for SIS 1 (the maximum possible extraction radii for these detectors), 
and « 6 arcmin for G1S2 and G1S3. For observation 3 we excluded the central 40arcsec to avoid the possibility of photon pileup. 
We used the sisclean and gisclean tools (with default values) to remove hot and flickering pixels. We filtered the data with the strict 
cleaning criteria outlined in the ASCA Data Analysis Handbook; however, we took the larger value of 7 GeV/c for the rigidity. We 
rebinned the spectral files so that each energy bin contained a minimum of 20 photons. We retained SIS data in the 0.5 to 9 ke V range 
and GIS data in the 1.5 to 9keV range. 

We accounted for the cross-calibration uncertainties of the three instruments by introducing a multiplicative constant for each 
detector in all of our fits. Relative to SISO, the SIS 1 detector normalization was always found to be within 2%, the G1S2 normalization 
was found to to be within 9%, and the G1S3 normalization was found to be within 15%. For any given observation, the normalization 
constants varied by ^ ±1% for different spectral fits. The background was measured from rectangular regions on the two edges of 
the chip farthest from the source (SIS data), or from annuli with inner radii > 6 arcmin (GIS data). These data were cleaned and 
filtered in the same manner as the source files. 

The resulting data files showed reasonable agreement between all four detectors. The most discrepant detector was SlSl, which 
also was the detector limited to the smallest extraction radius. This detector tended to show deviations from the other detectors for 
energies > 9 keV, and from the SISO detector for energies « 0.5-1 keV. The detectors were mostly in mutual agreement for the lowest 
flux observations. It is likely that the agreement could be further improved for observations located closer to the center of the chips 
(thereby allowing larger extraction radii), and if low galactic latitude dark sky observations in 1-CCD mode were available to use as 
background. 

RXTE DATA EXTRACTION 

In §1^ we present data from both pointed instruments on RXTE, the Proportional Counter Array, PCA, and the High Energy X-ray 
Timing Experiment, HEXTE. As we showed in §^ the large effective area of these instruments results in a data analysis approach 
that is dominated by the calibration uncertainty of these detectors (especially the PCA). In this appendix we summarize the major 
properties of both instruments and study their (inter-)calibration. All RXTE results obtained in this paper were obtained using the 
standard RXTE data analysis software, ftools version 4. 1 (including the RXTE -patch 4.1.1 and the correct accounting of the time- 
depe ndence of the PCA response; Jahoda, 1998, priv. comm.). Spectral modeling was done using XSPEC, version 10.00s (A rnaud 
1996). 

The PCA consists of five nearly identical co-aligned Xenon proportional counter u nits (PCUs) with a total effective area of about 



6500 cm^. The instrument is sensitive in the energy range from 2keV to ~ 60keV ( pahoda et al. 1996| ). We only used data where 
all five PC Us were turned on. Background subtraction was done in the same manner as for our RXTE Cyg X-1 observations (D ove 
et al. 1998). Specifically, a model using the rate of Very Large Events in the detector was used to estimate the background flux. 
The major uncertainty of this estimated background is due to activation of radioactivity in the detectors during SAA passages. Since 
this background component is present for about 30 minutes after the passage, we ignored data measured during these intervals. 
Furthermore, data were not accumulated at times of high electron contamination. The electron contamination is measured by a 
certain ratio of veto rates in the detectors, the so-called "electron ratio" (Jahoda et al., 1998, in preparation). As recommended by the 
RXTE Guest Observer's Facility (GOF), we excluded times during which the "electron ratio" was larger than 0.1 in at least one of 
the detectors. Note that the observed count rates from GX 339-4 are too high to allow the use of the newer background model made 
available by the RXTE GOF in 1998 June. 

For spectral fitting, we limited the energy range of the PCA data from 3 to 30keV and used version 2.2.1 of the PCA response 
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matrices. These matrices are newer than those used by us previously ( Dove et al. 1998 ), and they are primarily characterized by 
assuming a higher instrumental resolution (Jahoda et al., 1998, in preparation). Due to the large PCA count rate of GX 339-4 
(~ 800 cps) our observations are dominated by the remaining uncertainties in the detector calibration, and not by Poisson errors. 
Therefore, a good understanding of these uncertainties is necessary. 



Since the Crab spectrum is commonly assumed to be a featureless power-law (Poor & Seward 1974), at least over narrow energy 
ranges, the ratio between the fit to the Crab and the data can be used to deduce the systematic uncertainty associated with the detector 
caUbration. We therefore extracted a public domain spectrum of the Crab nebula and pulsar measured with the PCA on 1997 April 1. 
The Crab data were screened using the same criteria as those applied to our GX 339-4 data, except for that the "electron ratio" check 
was not applied since the background contributes only 0.6% to the total number of photons detected for the Crab. Modeling the 3 to 
30keV Crab data with an absorbed power-law resulted in a best-fit photon spectrum of the form 



13.3£"^ ''^^exp [-2.54 x lO^'cm-^bK^)] 



cm ^s 



keV" 



(Bl) 



where CTbf(£') is the en ergy dependent bound free absorption cross section for material of cosmic abundances as given by M orrison 
& McCammon (1983), and where E is the photon-energy measured in keV. For this fit to the Crab data, x^/dof = 168 /56. In Fig. ^ 
we display the ratio between the best-fit to the Crab and the data. Using this ratio plot, we deduced the systematic uncertainties of 
the detector (Table B9). Adding them in quadrature to the Poisson errors of our data, the of the above Crab fit was reduced to 
0.56. Note that the internal consistency of the detector calibrati on of the PCA app ears to be on the 1 % level, i.e., comparabl e 



Xred 



to that obtained for previously flown instruments like the Ginga LAC ( Turner et al. 1989[) or the ASC A GIS ( Makishima et al. 1 996 ), 



even though the much larger effective area of the PCA makes it necessary to include many "dirt effects" into the detector model. 

We caution that the statistical uncertainties of the fit parameter s presented in this work were derived u sing the above systematic un- 
certainties. It is questionable, therefore, whether the approach of Lampton, Margon & Bowyer (1976) to determine the uncertainties 
from the x^-contours can really be used since this approach makes use of the assumption of Poisson-type errors. The uncertainties 
given in this work should be taken with these caveats in mind, especially for those fits where the Xred values are very small (i.e., 
^ 0.5), and thus they should not be construed literally as '90% uncertainties'. 

The HEXTE consists of two cluster s of four Nal/Csl-phoswi ch scintillation counters, sensitive from 15 to 250 keV. A full de- 



scription of the instrument is given by Rothschild et al. (1998) Source-background rocking of the t wo clusters provides a direct 



measurement of the HEXTE background with measured long term syst ematic uncertainties of < 1% ([Rothschild et al. 1998 ). Al- 



though no other strong so urces in the field around GX 339-4 are known ( [Covault, Grindlay & Manandhar 1992 ; pouchet et al. 1993 
Trudolyubov et al. 1998), we extracted individual background spectra for both HEXTE cluster background positions to check for 



contamination of the spectrum from weak background sources. In all cases the background spectra differed by ^ 1 cps. Thus, we 
used the added background spectra from both cluster positions in our data analysis. We used the standard response matrices dated 
1997 March 20 and used data measured between 17 and 1 lOkeV. An analysis of the detector calibration similar to that performed for 
the PCA reveals that the HEXTE calibration is good on a level comparable to the PCA. Due to the much smaller effective area of the 
detector and due to the smaller flux from the source at higher energies, however, the HEXTE spectrum is completely dominated by 
the Poisson error of the data. Therefore we did not attempt to take the systematic calibration uncertainty into account. To improve 
the statistics of individual energy bins, we rebinned the raw (« 1 keV wide) energy channels by a factor of 2 for the energy range 
from 30 to 5 1 ke V, and by a factor of 3 above that. 



1.02 - 



1.00 - 



"5 0.98 - 



D 0.96- 



0.94 - 



0.92 



0.90 L 



10 

Channel Energy [keV] 
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Fig. B8. — Ratio between the power-law fit of eq. ( pH ) and the PCA Crab nebula data. The major deviations between the data and the model are 
due to the Xe L 1,2,3 edges at ~ 5. 1 keV, and due to the Xe K edge at 34.6 keV. 
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When modeling the spectrum of GX 339-4 from 3 to 110 keV as measured by the PCA and the HEXTE, the intercalibration 
betwe en the instruments is of some concern. Our experience with previous data and the response matrix described by D ove et al. 
(1998) indicated that the flux calibration of the HEXTE with respect to the PCA was off by about 25%, i.e., the derived HEXTE 
fluxes were ^75% of those found with the PCA. This deviation is mainly due to a slight misalignment of the HEXTE collimators 
which has not yet been taken into account in the HEXTE response matrix (Heindl 1998, priv. comm.). Using the new PCA response, 
the flux ratio now appears to be larger and was found to be 62 — 69% in our data. Extracting spectra with internal software used by 
the HEXTE hardware team produced spectra identical to those found using our extraction procedure. Therefore we do not believe 
this change in the flux calibration to be due to errors in the deadtime correction. To take this offset in the effective areas into account 
we modeled the spectra using a multiplicative constant which was set to unity for the PCA, and which was a fit parameter for both 
HEXTE clusters. Thus, all fluxes given below were measured with respect to the PCA. The maximum deviation of the HEXTE 
clusters relative to each other was found to be less than 8%. 

Apart from the flux calibration, however, even more crucial for our analysis is the question of how well the inferred spectral shapes 
agree for the two instruments. Our Crab fits show that the PCA results in a photon index of F = 2.187 (Eq. Bl), while our HEXTE 
fits gave r = 2.0 53. The generally accep ted value for the Crab photon index in the 1-100 keV range is F = 2. 10 ± 0.03. This value 
was adopted by Toor & Seward (1974)| in t heir analysi s of 28 different rocket flight measurements. There are indications that the 
spectrum softens to F ~ 2.5 above 150keV ( Jung 1989| ). Althoug h the absolute uncertainty of t he Crab flux in the 2-100 keV range 
has been estimated as large as 2%, and even larger below 2 keV ( N'0rgaard-Nielsen et al. 1994), the AF = 0.134 deviation between 
the PCA and the HEXTE best-fit Crab photon index is still very worrisome, and it is currently being studied by both the PCA and 
HEXTE instrument teams (Jahoda 1998, priv. comm.). 
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